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Molecular dynamics simulations have been used to study a collective atomic transport phenomenon 
by repeated Ar + irradiations in the Ti/Pt interfacial system. The ion-induced injection of surface 
atoms to the bulk, the ejection of bulk atoms to the top layers together with surface erosion is strongly 
enhanced by interfacial mixing. This process leads to a dense interfacial material, and broadening 
of the interface region. The process scales with the relative difference of the atomic masses. We 
find that surface roughening and interfacial mixing is strongly coupled via an enhanced counterfiow 
material transport normal to the surface which might be a novel surface damage mechanism. This 
cooperative phenomenon is active when the bilayer system is subjected to a high dose ion irradiation 
(multiple ion irradiations) and leads to surface cavity growth. 
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I. INTRODUCTION 

Low-energy ion bombardment (around some keV) is 
a widely used technique in surface analysis, sputter 
cleaning, depth-profiling, micromachining, TEM speci- 
men preparation etc. A tremendous experience on var- 
ious effects of ion irradiation on solid surfaces has been 
amassed [1]. There are also several theoretical descrip- 
tions of the process especially if we consider homogeneous 
material. The most visible effect of the ion bombard- 
ment, which can be studied easily experimentally (e.g. 
by means of TEM, STM, AFM) is the change in the to- 
pography of the surface. It is generally accepted that the 
surface morphology produced by subjecting the surface 
to a given ion fluence is the result of competition between 
various processes resulting in smoothing and roughening. 
The roughening can produce somewhat periodic struc- 
tures and/or rough structures which can be character- 
ized by various scaling lows (a detailed reference of the 
experimental results are given in ref. [2,3]). 

Molecular-dynamics (MD) simulations have also 
demonstrated the existence of the above mechanisms and 
proved to be capable to describe the development of vari- 
ous surface features during low energy ion sputtering and 
bombardment [4,6-10]. 

Three categories of atomic mass transport can be clas- 
sified when surfaces are subjected to ion beam irradiation 
[1]: (i) Atoms are ejected from the top layers, usually 
from the surface and srictly speaking this is what we call 
"sputtering" [2,3]. (ii) The displacement of atoms within 
their own layers creating vacancies, eroded surface, etc. 
(iii) Intermixing between layers [4] . 

These transport phenomena usually lead to the change 
of the surface morphology when the thermal spike [5] is 
located close to the free surface. It has been shown re- 



cently that in this special situation the surface damage 
mechanisms could be chategorized to four separate mech- 
anisms [6]: ballistic collisions, liquid How [4], microexplo- 
sions [11] and coherent displacements [12]. Recent exten- 
sive MD and experimental studies showed that the ion 
beam irradiation leads to surface modifications that are 
more extensive than just simple sputtering [6]. Ion irra- 
diation of interfacial metal bilayers are much less studied 
by MD [13-15] although the literature of interfaces in 
crystalline solids is rich [16]. 

In this paper we will report on MD simulation of ion 
irradiation induced alteration in the case of 4 Ti atomic 
layer/Pt substrate system. This system is specific in that 
sense that we place the interface close to the free surface 
(~ 10A) hence the thermal spike effects occur not only 
close to the free surface but also in the vicinity of the 
interface. The irradiation conditions were: ion energy 
1 keV at grazing angle of incidence (83°). Though the 
applied conditions can be considered gentle still strong 
mixing was observed. We will conclude that because 
of the presence of two unlike atoms enhanced material 
transport resulting in roughening and mixing occurs. We 
would like to show that in the case of a sample made of 
thin layers of different atoms with large difference in the 
atomic mass, not ion sputtering causes the surface erosion 
but it is rather induced by interfacial mixing. Morover 
the mixing of atoms with different atomic mass and the 
cavitation at the surface is strongly coupled via the ion 
induced atomic transport normal to the surface leading 
to a cooperative phenomenon. This process might be a 
novel surface damage mechanism never reported before. 

We irradiate a Ti/Pt system using several ion impacts 
consecutively. The reason for choosing such an irradia- 
tion mode is evident, since practically in all applications 
a series of ions hit the surface, even for low ion doses 



1 



[20]. We show through MD simulations of 1 keV Ar ions 
impacting the Ti surface of the Ti/Pt system that near 
surface interfacial mixing has a profound effect on the 
mass transport. 

II. SAMPLE PREPARATION AND SIMULATION 
TECHNIQUE 

Classical molecular dynamics simulations (MD) were 
used to simulate the ion-solid interaction using the PAR- 
CAS code developed by Nordlund et al. [17]. Here we 
only shortly summarize the most important aspects. Its 
detailed description is given in [17] and details specific to 
the current system is in a recent communication [15]. 

The sample consists of 37289 atoms for the interface 
(IF) system with 4 Ti top layers and a bulk which is 
Pt. The lattice constants for Pt is a « 3.92A and for 
Ti a « 2.95 and c « 4.68 A. At the interface (111) of 
the fee crystal is parallel to (0001) of the hep crystal and 
close packed directions are parallel. The interfacial sys- 
tem as a heterophase bicrystal and a composite object 
of two different crystals with different symmetry are cre- 
ated as follows: the hep Ti is put by hand on the (111) 
Pt bulk and various structures are probed and are put 
together randomly. Finally that one is selected which has 
the smallest misfit strain prior to the relaxation run. The 
remaining misfit is properly minimized below ~ 6% dur- 
ing the relaxation process so that the Ti and Pt layers 
keep their original crystal structure. The correspond- 
ing Ti-Ti and Pt-Pt interatomic distances are 2.89 and 
2.73 A at the interface. The Ti (hep) and Pt (fee) lay- 
ers at the interface initially are separated by 2.8 A and 
allowed freely to relax during the simulations. The varia- 
tion of the equilibrium Ti-Pt distance within a reasonable 
2.6 — 3.0 A interval in the interatomic potential does not 
affect the final results significantly. We find the average 
value of d s=s 2.65 A Ti-Pt distance in the various irra- 
diation steps and also in the nonirradiated system after 
a careful relaxation process. We believe that the system 
is properly relaxed and equilibrated before and between 
the irradiation steps. We found no vacancies at the in- 
terface or elsewhere in the system after the relaxation 
procedure. We visually checked and found no apparent 
screw dislocations, misorientations or any kind of unde- 
sired distorsions (in general stress-generator dislocations) 
at the interface or elsewhere in the system which could 
result in stressed induced rearrangements in the crystal 
during the simulations. The temperature scaled softly 
down towards zero at the outermost three atomic lay- 
ers during the cascade events [17]. The simulations are 
repeated in a much larger system consisting of ~ 72000 
atoms in order to rule out any kind of size effects. The 
simulation cell is of particular interest because it has been 
shown [15] that at the applied sputtering conditions (1 
keV Ar + with grazing angle of incidence) the thermal 
spike is typically formed at the interface and thus the 
emerging processes can be easily studied. 



The initial velocity direction of the impacting atom 
was 7 degrees with respect to the surface of the crystal. 
Further details are given in Ref. [15]. A series of up to 
10 ion impacts are employed in order to study the ef- 
fects of multiple ion irradiations, e.g. to simulate high 
dose effect. Between impacts the simulation time was 
<~ 20 ps. This time was found to be sufficient to reach 
a relaxed system after each ion impact (e.g. the tem- 
perature approaches zero). Thus direct interaction be- 
tween the cascade processes of various irradiation steps 
are excluded. It should also be remarked that the dam- 
age rate is extremely high due to the high dose effect and 
therefore some recovery process might occur between the 
irradiation steps. Although we did not find significant 
recrystallization effects up to several hundred ps. Ther- 
mally activated process might induce, however, further 
recovery in the system. We have done relaxation simu- 
lations on the irradiated sample at 300 K and found no 
recrystallization of the amorphised top layers in the ir- 
radiated zone up to several hundred ps. In any case, we 
do not think that the alloy phase would phase separate 
either at room temperature or on a longer time scale. We 
will demonstrate in the Discussion section that the the 
alloy phase is energetically competitive with the phase 
separated system (based on the cohesive energies). 

To obtain a representative statistics, the impact po- 
sition of the incoming ion is varied randomly within a 
10 x 10 A 2 area. We find considerable variation of mixing 
as a function of the impact position of the recoil. Using 
a 20 x 20 A irradiation area (with 80000 atoms or more) 
we checked that a further increase of the irradiated area 
does not affect the final results heavily, except that the 
larger the area, the larger the number of incoming ions 
needed to get the same extent of interfacial mixing (IFM) 
and surface roughening. This is because individual ions 
with impact positions far from each other do not result in 
the high dose effect (~ 10 15 Ar + ions/cm 2 ) we are inter- 
ested in in this article. Impact positions with sufficient 
proximity to each other (that is within the characteristic 
area ~ 5x5 A 2 , ~ 10 15 ion/cm 2 ) lead to the phenomenon 
called mixing induced enhancement of surface roughening 
(MIESR). The overlap of these individually eroded sur- 
face areas might lead to the macroscopic phenomenon of 
sputter removal when a sufficiently high dose of ions are 
bombarded to the surface in that case when sputtering 
is the only active mechanism for surface damage. How- 
ever, the final morphology of the sputtered surface might 
also be affected by MIESR if sufficiently large interfacial 
mixing occurs close to the free surface. 

We saw no further recovery process after 20 ps sim- 
ulations although we have extended several runs until 
500 ps. Thermally activated processes do not also play a 
significant role up to moderate temperatures (room tem- 
perature or so). Also we set the quenching rate to a value 
(5 K/Fs) which seems to be a reasonable choice. Reason- 
able variation of the quenching rate does not affect the 
cooling of the system and does not lead to the recovery 
of the original crystal order in the upper layers. 
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FIG. 1. The cross-sectional view of Ti/Pt after 5 Ar + 
irradiations (upper figure) at natural atomic mass ratio 
(niTi/rnpt ~ 0.25). The cross-sectional view of Ti/Pt af- 
ter 3 Ar + irradiations at the mass ratio of 0.06 (lower figure). 
The upper layer Ti atoms are shown with small circles. The 
z-axis corresponds to the depth and the x-axis to the horizon- 
tal positions in A . 



During the simulations we observe extensive depres- 
sion of the bombarded surface with respect to the non- 
sputtered surface. This depression did not result in a 
flat bottom, rather a rough surface had formed. We es- 
timate the extent of cavitation by counting the number 
of missing atoms in the cavity using the same procedure 
introduced for vacancies [15]. 

Vacancies were recognized in the simulations using a 
very simple analysis. A lattice site with an empty cell 
was considered to be a vacancy, where the diameter of 
the cell sphere is set to ~ 2A (the radius of a sphere 
around the relocated atom around its original position) 
which is about 60% of the average atomic distance in this 
system. We find this criterion is suitable for counting the 
number of vacancies. 

We used a many body potential, the type of an 
embedded-atom- method given by Cleri and Rosato [18], 
to describe interatomic interactions. This type of a po- 
tential gives a very good description of lattice vacancies, 
including migration properties and a reasonable descrip- 
tion of solid surfaces and melting [18]. Since the present 
work is mostly associated with the elastic properties, 
melting behaviors, surface, interface and migration en- 
ergies, we believe the model used should be suitable for 
this study. The Ti-Pt interatomic potential of the Cleri- 
Rosato [18] type is fitted to the experimental heat of mix- 
ing of the corresponding alloy system. Further details are 
also given elsewhere [15]. 

An atom is labeled mixed if it moved beyond the in- 
terface by more then d/2 w 1.4A along the z-axis, where 
d is the interfacial Ti-Pt distance. The broadening at 
the interface (or interfacial mixing, IFM) is calculated as 
follows: 
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FIG. 2. Cavitation (cratering) in the Ti/Pt system after 
the 5th bombardment. The xy-plot of the Ti/Pt system seen 
from above (from the free surface). 



where N m i x is the number of mixed atoms and Azi is 
the distance of the mixed atom i from the interface (z- 
axis depth position), b is a useful quantity for compar- 
ison with measurements only in that case when Azi is 
larger then the interlayer distance. The reason is that 
the measured broadening is dominated by those inter- 
mixed atoms which move beyond the interface by more 
than the inter-layer distance (2nd layer mixed atoms). 
Using this condition we get a reasonable agreement with 
Auger depth profiling [20]. 



III. RESULTS AND DISCUSSION 

The detailed analysis of the mixing process in the 
Ti/Pt interfacial system after a single ion irradiation pro- 
cess has been given recently [15] . Here we restrict oursclfs 
to the study of the effect of multiple ion impacts on mix- 
ing as well as on surface erosion. 
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FIG. 3. The Cross-sectional view of pure hep Ti (upper 
figure) and fee Pt (lower) after 5 Ar + irrad iations. 



In FIGs (1)- (2) the xz cross-sectional and xy top 
plots of the Ti/Pt system are shown after the 5th irradi- 
ation. It is evident that a considerable amount of mate- 
rial has been removed from the irradiated area. Since the 
height of the non-irradiated surface did not change apart 
from some adatoms the material removal produces a fea- 
ture like a crater. It should be added, however, that the 
feature formed has a rough surface. The volume of the 
crater is much larger than which is expected considering 
only sputtered atoms. The number of sputtered atoms 
is less by an order of magnitude then the number of in- 
termixed Ti atoms which are pushed into the Pt bulk. 
These atoms cause the considerable interface broadening 
visible in FIGs (1) and (4). Indeed, the majority of 
removed atoms are the intermixed Ti atoms (injected to 
the bulk). In FIGs (3) we show the corresponding pure 
components after the 5th bombardment and see no ap- 
parent surface morphology development. This also shows 
that the degree of the material removal by sputtering is 
much weaker than the amount of the missing atoms in 
the case of the multilaycrcd material. Since the differ- 
ence between the pure material and the Ti/Pt system is 
that in the case of the latter a two component system 
forms due to ion sputtering, we should conclude that the 
presence of the second component seriously affects the 
mass transport and thus most likely leads to morphol- 
ogy development as well. Let us check more closely the 
correspondence between IFM and mass transport (which 
is measured by the volume of the apparent crater). Our 
primary intention is then to point out that the enhance- 
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FIG. 4. Upper figure: broadening at the interface (A) as a 
function of the irradiation steps. Lower figure: the number 
of removed (missing) atoms in the surface cavity as a func- 
tion of the irradiation steps. Note that most of the removed 
atoms are injected to the bulk and only a small portion of 
them (~ 10%) are sputtered to the vacuum. The number of 
vacancies in the whole sample are also given. 



ment of surface cavitation in Ti/Pt is due to interfacial 
mixing. MIESR is strongly connected to a two-way ma- 
terial transport (with countcrflow motion) which might 
be understood as a radiation-induced flow of Pt atoms to 
the surface and the injection of Ti to the bulk. 

We see the strong enhancement of IFM due to the re- 
peated irradiation of the surface in FIG (4) where the 
broadening at the IF is plotted against the number of 
irradiation steps (upper figure). Anomalously large in- 
termixing in metal bilayers have recently been seen as 
well by Buchanan et al. [21]. Surprisingly, we see the 
saturation of mixing (broadening) at the 5th step. The 
magnitude of mixing at saturation shows no dependence 
on heat of mixing as it was obtained in the first irradia- 
tion step [15]. Hence MIESR is not governed by thermal 
spikes. 

In the lower part of FIG (4) the number of "miss- 
ing" atoms in the top layer (vacancies) is plotted against 
the number of irradiation steps. The number of miss- 
ing atoms is proportional roughly to the volume of the 
surface cavity which are mostly injected to the bulk. In- 
terestingly these features coincide with each other: the 
number of intermixed Ti atoms (the "missing" atoms) 
saturates also at the 5th step. 

The volume of the cavity rapidly increases in the first 
three irradiation steps and a crater similar to that of re- 
ported in earlier studies by Bringa et al. at higher impact 
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FIG. 5. The number of removed atoms in the top layers 
(the extent of cavitation) as a function of the atomic mass ra- 
tio (5 = mtop/mpt) in the third irradiation step. The straight 
line is a linear fit to the data points, and the error bars denote 
standard deviations. 



energies with single ion bombardment, forms [8]. In these 
cases cratering is attributed to the high deposited energy 
close to the free surface at certain nonchanneling direc- 
tions. In our case, however, after repeated irradiation of 
the surface, there seems to be a strong correlation be- 
tween interfacial intermixing and material removal from 
the surface (material injection to the bulk). Pt atoms 
intermixed to the top layers, however, do not occupy the 
empty sites in the cavity. Rather they form a dense phase 
at the surface of the cavity. The concentration of atoms 
in the irradiated region (10 x 10 x 10 A 3 ) increases from 
the initial 6.1 x 10 22 atom/cm 3 to 6.8 x 10 22 atom/cm 3 
in the 5th irradiation step. The atomic concentration 
n in TiPt alloy is ~ 7.2 x 10 22 atom/cm 3 , therefore the 
increase of n in the irradiated region indicates the forma- 
tion of a broad alloy phase at the interface. However, the 
saturation of the curves in FIGs (4) indicates that the 
simultaneouss process of intermixing and surface cavi- 
tation is limited. A possible control factor can be the 
increase in the number of the undercoordinated surface 
atoms which is always unfavored. Another factor could 
be the saturation of the broadened interface. No more 
intermixing of Ti atoms is observed possibly due to the 
denseness of the alloy phase formed in the original Pt 
bulk which does not permit the further transport of Ti 
atoms to the Pt bulk. Anyhow the maximal broadening 
of the IF seems to be limited to b s=a 10A. This restric- 
tion on broadening seems not to be affected by chemical 
forces since it was found that the heat of mixing did not 
affect the IFM [15]. 

In order to understand more fundamentally the mech- 
anism of MIESR we vary the atomic mass of the top layer 
constituent (m top ), keeping the same potential which can 
be done easily in MD simulations [17]. According to our 
expectations we find the suppression of IFM and surface 
cavitation when m t0 p/mpt — > 1. 

We plot in FIG (5) the number of Ti atoms which 
was removed from the top layers by MIESR (cavitation) 
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FIG. 6. The cohesive energy (eV/atom) of Pt and Ti atoms 
as a function of their depth position (A ) in the 10th irradia- 
tion step. The dotted line corresponds to the depth positions 
of the interface and the free surface is at depth position 0. 
The cohesive energy in pure Pt and in Ti is ~ 5.8 eV/atom 
and ~ 4.8 eV/atom, respectively, which should be compared 
with the value of 6.3 eV/atom in TiPt alloy. 



as a function of the atomic mass ratio obtained from 
simulations in the third irradiation step. In this step we 
obtained extensive mass transport especially for lighter 
atoms (mtop ~ 2 — 11 g/mol, S « 0.01 — 0.05, see also FIG 
(1)). Notable feature is in lower FIG (1) the enhanced 
cavitation with small mass ratio in the third irradiation 
step. 

It can obviously be seen that the lighter the atoms of 
the top layers, the larger the extent of the mass transport. 
At the same time we see no obvious effect of heat of 
mixing on IFM [15]. It is therefore quite probable that 
IFM and surface roughening is governed primarily by the 
atomic mass difference between the constituent atoms of 
various layers. Mass transport from the surface layer 
is strongly enhanced for multilayers with light atoms in 
the thin top layer. This finding supports our previous 
conclusion [15] that the ballistic mechanism governs the 
IFM. 

This is interesting because we are not aware of similar 
results where the mass transport is strongly coupled to 
interfacial mixing. Moreover the bulk process (mixing) 
and surface roughening arc coupled via the enhanced ma- 
terial transport normal to the surface leading to a coop- 
erative phenomenon. The importance of this finding is 
to point out that the surface of materials which include 
more components close to the surface might behave in a 
different way than the pure components. The surface ero- 
sion is governed not primarily by sputtering but rather 
by mixing. A specific feature of this system is that we 
found no effect of chemical forces (heat of mixing) on 
IFM [15]. The process is governed by the large difference 
of atomic mass between Ti and Pt. 

In order to understand the energy saldo of mixing we 
plot in FIG (6) the cohesive energy (binding energy) 
of various atoms as a function of their depth position. 
The data points are collected from the 10th irradiation 
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step. The atomic cohesive energy E co h is averaged with 
that of every atom's nearest neighborsrs to reduce the 
atom-level scatter in the values. It can obviously be seen 
that the cohesive energy of Pt atoms decreases when the 
atoms move towards the surface while Ti atoms gain en- 
ergy when mixed to the Pt bulk. The energy gain per Ti 
atom is ~ 1 eV on average while approximately the same 
energy loss occurs for the intermixed Pt atoms. Above 
the depth position of ~ —30 A the widening of the co- 
hesive energy "window" can also be seen for Pt. This 
is due to the formation of an alloy phase at the broad- 
ened interface. Close to the free surface no features of 
alloy phase formation can be seen, e.g. the E co h values 
are rather close to the atomic value of pure Ti. Also the 
E co h of the intermixed Pt atoms decreases at the sur- 
face which is due to their decreased coordination number 
(< 12) with respect to the bulk value of 12. 

Finally it should also be mentioned that the ion in- 
duced surface cavity growth might have many common 
features with the growth of cavities found on boundaries 
in stressed polycrystals or with cracking and in general 
with various fracture processes at heterophase interfaces 
[16]. This is because we do not rule out the possibility 
of a cavity growth mechanism induced by a stress field 
generated by the thermal spike. This is, however, beyond 
the scope of the present article and needs further studies. 

IV. CONCLUSIONS 

We have shown that the mass transport for the stud- 
ied double layer system is different from that of the cor- 
responding pure elements. We have also shown that the 
observed apparent cavitation is in correlation with the 
interfacial mixing and it increases with increasing mass 
difference of the constituents. The most important conse- 
quence of this finding is that a reliable theory describing 
the surface morphology development of a system contain- 
ing two components should also account for the mixing 
near the surface. Presently available theories are miss- 
ing this point. This novel surface damage mechanism 
might be applicable for those multicomponent systems 
(e.g. thin films) where the material is subjected to low 
energy ion-irradiation and the interface is close to the 
free surface. The mechanism is regressive ly reliable for 
systems with mass ratio approaching 1. 
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